Strong Upper Limits on Sterile Neutrino Warm Dark Matter 
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Sterile neutrinos are attractive dark matter candidates. Their parameter space of mass and 
mixing angle has not yet been fully tested despite intensive efforts that exploit their gravitational 
clustering properties and radiative decays. We use the limits on gamma-ray line emission from the 
Galactic Center region obtained with the SPI spectrometer on the INTEGRAL satellite to set new 
constraints, which improve on the earlier bounds on mixing by more than two orders of magnitude, 
and thus strongly restrict a wide and interesting range of models. 

PACS numbers: 95.35.+d, 13.35.Hb, 14.60.St, 14.60.Pq 



Introduction. — The existence of dark matter is cer- 
tain, but the properties of the dark matter particles 
are only poorly constrained, with several attractive but 
rather different candidates. One of these, sterile neutri- 
nos, would be a very plausible addition to the Standard 
Model [H, 0, 13- If their masses were in the range 
~ 0.1 — 100 keV, they would also act as "warm" dark 
matter [1, H, LI , which could be in better accord with ob- 
servations than standard "cold" dark matter candidates. 
Even if sterile neutrinos are not a dominant component 
of the dark matters they may still exist and cause other 
interesting effects [E@], such as pulsar kicks 0], and may 
affect reionization [8| . It is therefore important to deeply 
probe the sterile neutrino parameter space, as defined 
by the mass m s and mixing sin 2 26 with ordinary active 
neutrinos, and shown in Fig. [T] 

One means of testing sterile neutrino dark matter mod- 
els is through cosmological searches, which rely on the 
effects of sterile neutrino dark matter on the large-scale 
structure of gravitationally-collapsed objects. While re- 
cent results based on the clustering of the Lyman-a forest 
and on other data have been interpreted as lower limits 
on the sterile neutrino mass of up to about 10-13 keV, in- 
dependent of the mixing angle [9(] , these constraints may 
be weakened depending on the sterile neutrino produc- 
tion model (e.g., Ref. [l0(). 

Another means of constraining sterile neutrino dark 
matter is through their radiative decay to active neutri- 
nos, v s — ► v a + 7. These decays produce mono-energetic 
photons with E~, = m s /2. While the decay rate is exceed- 
ingly slow due to the tiny active-sterile mixing, modern 
satellite experiments can detect even these very small 
x-ray/gamma-ray fluxes, and such a signal could specif- 
ically identify a sterile neutrino dark matter candidate. 
The signal from nearby dark matter halos is line emission 
and the cosmic signal from all distant halos is broad- 
ened in energy by the integration over redshift. There 



are limits obtained using the Cosmic X-ray Background 
(CXB) data @, [ll| and, at lower masses, stronger lim- 
its using data from a variety of nearby sources (see e.g., 
Refs. [H, Qj3] and references therein). 

It is important to improve on both the cosmological 
and radiative decay constraints; despite intensive efforts, 
viable models that match the observed dark matter den- 
sity still remain. In fact, it has recently been empha- 
sized [1, |i| that some models may extend to regions of 
the parameter space far from the earliest and simplest 
models [3] to much smaller mixing angles. We calculate 
the gamma-ray flux from dark matter decays around the 
Milky Way center and compare this to the limits on the 
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FIG. 1: The sterile neutrino dark matter mass m 3 and mixing 
sin 2 26 parameter space, with shaded regions excluded. The 
strongest radiative decay bounds are shown, labeled as Milky 
Way (this pap er), CXB II ill, and X-ray Limits (summarized 
using Ref. [12J ; the others [13J are comparable) . The strongest 
cosmological bounds 0] are shown by the horizontal band 
(see caveats in the text). The excluded Dodelson-Widrow 0] 
model is shown by the solid line; rightward, the dark matter 
density is too high (stripes) . The dotted lines are models from 
Ref. [l4|], now truncated by our constraints. 
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line emission flux from the INTEGRAL satellite. The 
high sensitivity and spectral resolution of the available 
data enable us to derive new and very stringent con- 
straints. For masses above 40 keV, this improves on the 
CXB constraints ll| on the mixing angle by more than 
two orders of magnitude. 

INTEGRAL Gamma-Ray Line Search. — Tee- 
garden and Watanabe have reported results from a 
search for gamma-ray line emission from point and dif- 
fuse sources in the Milky Way [3], using the SPI spec- 
trometer on the INTEGRAL satellite [16(. In the en- 
ergy range 20-8000 keV, they tested for lines of intrin- 
sic width 0, 10, 100, and 1000 keV. The additional line 
width due to instrumental resolution increases over the 
above energy range from ~ 2 to 8 keV full-width half- 
maximum (FWHM). As expected, their analysis recov- 
ered the known astrophysical diffuse line fluxes at 511 (l7| 
and 1809 keV [l8j], and no others, validating their proce- 



dures [15j. The principal advantages of the SPI instru- 



ment for a sterile neutrino decay search are its wide field 
of view and excellent energy resolution. For sterile neu- 
trino decays in the Milky Way halo, the line width due to 
virial motion is ~ 10 -3 , which is therefore small enough 
to be neglected. 

Two large-scale regions around the Galactic Center 
were considered, with angular radii of 13° and 30°, and 
exposures of 1.9 and 3.6 Ms, respectively. The 24° col- 
limated field of view was used without the coded mask 
image reconstruction and the corresponding limits on the 
flux from an unknown line emission were derived by de- 
convolving an assumed sky brightness distribution (ei- 
ther a Gaussian with 10° FWHM for the former or flat 
for the latter region) and the wide angular response of 
the collimator. To improve the sensitivity to line emis- 
sion specifically from these regions, the average flux away 
from the Galactic Center Region (angular radii of > 30°) 
was subtracted from the flux from inside the Galactic 
Center Region. This procedure cancels almost all of the 
instrumental backgrounds. This also cancels all of the 
cosmic signal and part of the halo signal, and a careful 
calculation of the latter effect is taken into account in our 
analysis. For the Galactic Center Region, the 3.5er limits 
on narrow line emission are < 10 -4 photons cm -2 s _1 for 
the full range of energies. The actual energy dependence 
of the limiting flux, !Fu m (E), is more complicated, and 
we took this into account (leading to the slightly jagged 
edge of our exclusion region) . 

Milky Way Dark Matter Decay Flux. — To turn 
the INTEGRAL limits on generic line emission into con- 
straints on sterile neutrino dark matter, we calculated 
the expected gamma-ray emission from the decay of ster- 
ile neutrinos in the Milky Way (the INTEGRAL limits 
also strongly constrain certain decays of GeV-mass dark 
matter models flij]). For a long-lived decaying sterile 



neutrino with lifetime r and mass density p — m s n, the 
intensity [201 ] (number flux per solid angle) of the decay 
photons coming from an angle ?/> relative to the Galactic 
Center direction is 



47rm s r 

where the dimensionless line of sight integral, 
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is normalized at the solar circle, with R sc — 8.5 kpc and 
p sc = 0.3 GeV cm -3 (these cancel later). While t ma x 
depends on the adopted size of the halo, contributions 
beyond the scale radius of the density profile, typically 
about 20-30 kpc, are negligible. 

The sterile neutrino radiative lifetime r is 
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where we have chosen the Dirac neutrino decay life- 
time 21|; for the Majorana case, which may be favored, 



the lifetime is 2 times shorter, which would lead to more 
restrictive constraints. The prefactor in Eq. ([T]) can then 
be expressed in terms of the mass and mixing of the ster- 
ile neutrino, 
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The number flux of photons at energy E y = m s /2 is 
obtained by integrating the intensity, Eq. |T]), over the 
field of view, 



A<> 



PscRsc 

47rm,T 



dfl Jty) , (5) 



where the solid angle is A57 = 27r(l — cos ip). 

The dark matter distribution of the Milky Way is not 
perfectly known [i^] , though the variations between mod- 
els make little difference for dark matter decay, since the 
density appears only linearly in the calculations (unlike 
for dark matter annihilation, where it appears quadrat- 
ically). A trivial lower bound for the integral in Eq. §5§ 
can be obtained by taking the dark matter density to be 
constant within some radius from the Galactic Center, 
which we take to be R sc . Then the line of sight and field 
of view integrals are just multiplications: using Eq. 
the former is ~ 2, and since Afl ~ 0.16 for ip = 13°, the 
latter is f AQ dfl J(ip) ~ 0.3. 

For realistic dark matter density profiles, the field of 
view integral in Eq. ([5]) will be larger, since the density is 
larger (though more uncertain) in the central region. We 
calculated this for the Navarro-Frenk- White (NFW) [13|, 
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Moore 24|, and Kravtsov [2 5f p rofiles, which are all com- 
monly used (see also Ref. |26(). These are normalized 
with p(R sc ) = 0.30, 0.27, and 0.37 GeV cm" 3 , respec- 
tively. These slight differences in normalization compen- 
sate the different slopes at inner radii so that the masses 
enclosed at outer radii are the same [i^] • In the left panel 
of Fig. [2j the thin lines show Jfy) as a function of the 
angle ip for each profile; in the right panel, the corre- 
sponding thin lines show these integrated over the field 
of view (up to the angle ip), as in Eq. ([5]). These results 
take into account the variation of density with position, 
and also the contribution from halo dark matter beyond 
the solar circle on the other side of the Milky Way. Note 
that all three profiles have similar values of f. Q dtt J{ip), 
since the large field of view de-emphasizes the inner radii 
where the differences between the profiles are the largest. 

Constraints on Sterile Neutrinos. — As noted 
above, the INTEGRAL limits on line emission from the 
Galactic Center region are obtained by subtracting the 
average flux outside this region (ip > 30°) from the flux 
inside this region (ip < 13°), which must be taken into 
account in our analysis. To be conservative, we consid- 
ered the maximum effect of this subtraction by fixing the 
intensity outside the Galactic Center region to its value 
at ip = 30°. (In fact, it is smaller at larger angles.) In 
terms of our equations, this is 



d£l [J(iP) - J(30°)] 



(6) 



In the right panel of Fig. [2] our results for the integrated 
J{4>) — 17(30°) are shown by the thick lines. The effect of 
this subtraction correction is not large, less than a factor 
of 3 at ip = 13° for all three profiles. In addition, the 
INTEGRAL flux limits of Ref. [H| for an angular region 
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FIG. 2: Left panel: The line of sight integral J(ip) as a 
function of the pointing angle tp with respect to the Galac- 
tic Center direction for the three different profiles considered 
(Kravtsov, NFW, and Moore, in order of solid, dashed and 
dotted lines) . Right panel: Integrals up to the angle ip of J7"(V0 
(thin upper lines) and J(ip) — J7(30°) (thick lower lines). The 
grey line at 13° marks the field of view for the INTEGRAL 
flux limit, and we chose J dQ, [^T(tp) — JT(30°)] — 0.5 as a 
conservative value for our subsequent constraints. 



of ip < 13° assume that the line emission intensity follows 
a two-dimensional Gaussian with FWHM of 10°, while a 
flat-source profile would yield somewhat weaker limits. 

To shield our results from such uncertainties associated 
with the distribution of dark matter in the Milky Way, 
including whether warm dark matter profiles are less cen- 
trally concentrated than cold dark matter profiles, we use 
a rather conservative value, f. n dfl [J(ip) — ^(30°)] ~ 
0.5, in our subsequent calculations. Our results can be 
easily rescaled for a different value and our limits should 
improve as the amount of data increases in time. 

While we have presented our results for the region 
within 13° of the Galactic Center, there are also flux 
limits for an angular region of ip < 30° and an assump- 
tion that the intensity is constant in angle [15j . The flux 
limits for ip < 30° are ~ 3 times weaker than those for 
ip < 13° [15j. However, as shown in the right panel of 
Fig. [21 the sterile neutrino decay flux, which is propor- 
tional to J An dn [J{ip) - J(30°)], is ~ 2-3 times larger 
for ip < 30° than for ip < 13°, compensating the lower 
sensitivity. Thus our results are rather robust against the 
choice of angular region used and other assumptions for 
analyzing the INTEGRAL limits. 

With these detailed results on the sterile neutrino dark 
matter distribution, we define constraints in the param- 
eter space of mass and mixing. The expected line flux at 
E-y = m s /2 from dark matter decay, which depends on 
m s and sin 2 29, should not exceed the INTEGRAL limits 
(for 3.5ct), i.e., Tu m > AT S , or 
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Substituting Eq. and J AQ dtt [J{ip) - J(30°)} ~ 0.5 
yields our result in Fig. [I] The boundary of the excluded 
region is jagged on the left due to the actual energy 
dependence of the limiting flux, Tu m (E) (see Fig. 9 of 
Ref. 15j). The energy range available with the SPI in- 
strument causes the sharp cut-off at m s — 40 keV. Our 
constraint is coincidentally in line with prior constraints 
at lower masses using the x-ray emission from nearby 
sources. There is only a narrow gap, m s ~ 20-40 keV, 
in which the best available mixing constraints are sub- 
stantially weaker. The constraints shown in Fig. [1] as- 
sume that sterile neutrinos comprise all of the required 
present-day dark matter, but the limits at large mass are 
so stringent that they would provide strong limits even 
on sterile neutrinos that were only a fraction of the dark 
matter. 

Conclusions. — Sterile neutrinos require only a min- 
imal and plausible extension of the Standard Model [l], 
& & 0] and can solve problems in reconciling the ob- 
servations and predictions of large-scale structure 0, 0] . 
Despite intensive efforts on setting constraints, there are 
still viable sterile neutrino dark matter models over a 
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wide range of mass m s and mixing sin 2 28; the focus 
is now at larger mass and smaller mixing than consid- 
ered in the earliest and simplest models 0]. In this 
region, the models are very challenging to test, either 
through their differences in clustering with respect to 
cold dark matter candidates or their astrophysical ef- 
fects III . 

mi 




.or through their very small radiative decay 
H, lU or laboratory tests @ . 



rates 

Teegarden and Watanabe [15| presented the results of 
a sensitive search for line emission in the Galactic Cen- 
ter Region, using data from the SPI spectrometer on the 
INTEGRAL satellite [16[ . Based on a simple and conser- 
vative calculation of the expected gamma-ray flux from 
sterile neutrino dark matter decays, we have used these 
limits to set new and very strong constraints on sterile 
neutrino parameters, as shown in Fig. [1] The large-mass 
region is now very strongly excluded, improving on the 
previous CXB mixing constraints llj by more than two 
orders of magnitude. At fixed m s , the boundary in sin 2 28 
is defined by the 3.5a exclusion; using Eqs. ||7J) and ((4}, 
it is easy to see that points with sin 2 29 values ten times 
larger than at the boundary are excluded by a nominal 
35ct, and so on. On the scale of the figure, any reasonable 
further degradations in the conservatively-chosen inputs 
would not be visible. We anticipate that it will be possi- 
ble to extend our constraints, in particular going to lower 
masses, by dedicated analyses of the INTEGRAL data, 
which we strongly encourage. If the sensitivity of this 
and other techniques can be improved upon, then it may 
be possible to definitively test sterile neutrinos as a dark 
matter candidate. 
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